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Abstract 

We study the coherent final state interaction of an energetic parton produced in AA 
collisions caused by the change in the cutoff scale and running coupling constant from 
the vacuum to QGP. We demonstrate that the contribution of this new mechanism to the 
energy loss may be of the order of magnitude of the induced gluon radiation. However, 
an accurate evaluation of this medium effect is a difficult task because there is a strong 
cancellation between the cutoff and running coupling constant effects. The uncertainties 
in the contribution of the coherent final state interaction restrict strongly the accuracy of 
jet tomographic analyses of the matter density produced in AA reactions. 



1. Introduction. In recent years there has been much work done on the energy 
loss of fast partons in hot QCD medium due to gluon radiation induced by multiple 
scattering (for a review, see j^). This is of great importance for understanding final state 
interaction in hard reactions in high energy nucleus-nucleus collisions which are under 
active investigation at RHIC, and will be studied in future experiments at LHC 

The theoretical calculations show that the energy loss in quark-gluon plasma (QGP) 
exceeds considerably the one in hadronic medium P, ||, ||. Since gluon radiation softens 
the parton fragmentation functions of energetic partons produced in hard reactions in 
the initial stage of AA collisions it should lead to significant suppression of the high pt 
hadronic spectra in AA collisions with respect to pp collisions (so-called jet quenching) 
if a hot QGP is formed |^. Such a suppression was indeed recently discovered by the 
PHENIX experiment [0 at RHIC for tt^ spectra atp^ < 4 GeV in central Au+Au collisions 
at y/s = 130 GeV. Because the energy loss is sensitive to the density of hot medium it looks 
quite natural to use experimental data on high p^ spectra for jet tomographic analysis of 
the matter density produced in AA reactions p, ^ |10 . 



To understand the range of uncertainty in the jet tomographic analyses it is important 
to study other possible final state interaction effects in jet production. One mechanism 
of potential interest is the in-medium modification of the parton cascade without gluon 
exchanges between the fast partons and thermal partons (this mechanism we call the 
coherent final state interaction (CFSI)). The reason is evident: jet splitting in vacuum is 
the major mechanism of the energy loss of energetic partons, and if the medium affects the 
parton cascading one can expect a significant modification of the fragmentation functions. 
For example, such a modification should inevitably arise as a mass effect due to different 
infrared cutoff scales in the vacuum and QGP. Another obvious source of the CFSI is 
the in-medium modification of the running coupling constant as(fc).Q Although at large 
virtualities the running coupling constants in the vacuum and QGP are close to each other 
this should not be the case at low k where different background environments in which the 
gluon bremsstrahlung occurs can lead to a difference in the running coupling constants in 
these two cases. Note that both these medium effects should lead to transition radiation 
very similar to that of photon radiation in QED. The purpose of the present work is to 
address the CFSI for RHIC conditions within a simple model for gluon radiation which 
will be discussed in detail below. 

2. Cutoff scales and running coupling constants. Let us first discuss the 
magnitudes of the cutoff scales for parton splitting in the vacuum and QGP. In the QCD 
vacuum the natural cutoff is the inverse gluon correlation radius ~ 0.8 — 1 GeV 



111 , 0. Namely at such a virtuality scale the perturbative cascade stops in the Monte 
Carlo programs like JETSET, and the string fragmentation takes over. Note also that 
introduction of the effective gluon mass mg^^ ~ (hereafter we use the index v for 
vacuum quantities, and for the plasma quantities below we use the index p) allows one 
to describe the HERA data on the low-a; proton structure function [Q. The analysis of 
inclusive radiative decays of the J/\I' and T [|I^ also gives „ ~ 0.7 — 1.2 GeV. 

In the QGP phase the nonperturbative fluctuations are suppressed, and the natural 



^Note that these CFSI effects differ from the coherent double gluon exchanges which are usually 
included in the induced gluon radiation to insure the unitarity . 
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cutoff for radiation of transverse gluons, which propagate through QGP as quasiparticles, 
is the thermal gluon mass. At high temperature it reads m^ ^ = ("^ + "H^) • 
analysis of the results of the lattice calculations show that in the temperature range T ~ 
(1 — 3)Tc (Tc ~ 170 MeV is the temperature of the confinement phase transition), which is 
of relevance to AA collisions at RHIC, the nonberturbative effects are still important ||15|| . 
Using a quasiparticle picture with massive gluons and quarks in the above temperature 
window the authors of Ref. [l^ obtained nigp ^ 0.4 GeV, and rrig^p ~ 0.3 GeV. Thus we 
see that there is a considerable difference in the cutoffs in the vacuum and QGP. 

Let us recall now the situation with the running coupling constant at low k. There 
are some indications that in the vacuum nonperturbative effects stop the growth of the 
running coupling constant at A; < /c^ ~ 1 GeV [jl6|, |T^, |18|, |T^. Phenomenologically the 
magnitude of ag^v sX k <kc can be estimated from, say, the analysis of the heavy quark 



energy losses which gives [17 



2 GeV (y (k) 

rfA;^^^ 0.36 GeV. (1) 
TT ^ ^ 

For the simplest prescription with frozen a^^t, a.t k < k^ (the so-called F-model [|17|)5 using 
the one-loop expression at k > kc one can obtain from (|l]) as,v{k < k^ = a{.^'^ ~ 0.7, and 
fee ^ 0.82 GeV. These values are for Aqcd = 0.3 GeV. 

Unfortunately, at present there is no accurate information on as{k) for gluon emission 
from a fast parton in QGP. Available pQCD calculations are performed in the static 
limit (see for example pO|, |2l|, ^ and references therein). The running coupling constant 
obtained in ^ has a pole at k/ Kqcd ~ 3 at T ~ 250 MeV. Thus, in pQCD, even 
for the static case, the situation with fc-dependence of the in-medium running coupling 
constant at low k is unclear. On the other hand, analysis of the lattice results within 
quasiparticle model gives the thermal with a smooth T-dependence, and at T ~ 250 
MeV as ~ 0.5 fl^. In the present paper in the absence of accurate information on the 
in-medium running coupling constant for fast partons we perform calculations using the 
above F-model with different values of al'^„. 

3. Evaluation of coherent medium correction to gluon spectrum. Let us now 
discuss technical aspects of our analysis of the CFSI. We consider the gluon radiation 
from a fast quark (the generalization to the radiation from a gluon is trivial). We neglect 
multiple emission and consider only the leading order splitting q — >■ gq. It is reasonable 
since the effect is dominated by the gluons with small transverse momenta A; < 1 — 2 GeV. 
Note that we choose the z-axis along the momentum of the initial fast parton, so for the 
central rapidity region in AA collisions our L is the ordinary transverse distance between 
the jet production point and the boundary of QGP. 

We consider a fast quark with energy Ei produced at ^ = which eventually splits at 
some z > into the gluon and final quark with the energies Eg = xEi and Ef = {l — x)Ei 
respectively. The corresponding matrix element can be written in the form (below for 
simplicity we drop color factors) 

/■oo r _ 

T = ij^ dzj dpgijfip, z)YA^{p, z)ij,{p, z) , (2) 

where ipijlz, p) are the wave functions of the initial and final quarks, and is the wave 
function of the emitted gluon, p is the transverse coordinate. In (0) we do not explicitly 
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indicate the z and k dependence of the running couphng constant g. We evaluate the 
matrix element for small emission angles. Then, at high energies Ej ^ the quark 
wave functions using the ordinary light-cone spinor basis can be written as 

^Ap^ 2) = exp{iEjz)Uj(j)j{z, p) , (3) 

where the operator tjj reads 



, , ap + dm,, , , , 

|V2i^,+ ^ \x,. (4) 

Here Xj is the quark spinor (normalized to unity), cx = 7^7, (3 = 7°, and p = —iV±. The 
gluon wave function can be written in a form similar to (|^) (up to an obvious change of 
the spin operator). The transverse quark wave function (j)j{p,z) entering (^ is governed 
by the two-dimensional Schrodinger equation in which z plays the role of time 

»^^^.(..). (5) 

A similar equation holds for the gluon wave function. 

Without a loss of generality we can take for the initial quark the plane wave state in 
the p-plane and set Pj = 0. Then all the transverse wave functions can be written as 



[z, p) = exp 



p? + ^?(0) 



2E, 



(6) 



Eventually, the p-integration in (^ will give Pg + Pf = Pj = 0. Note that since the 
quark mass is of only marginal significance (for light quarks) in the gluon radiation we 
will neglect the z-dependence of ruq and use the same quark mass in the vacuum and 
in QGP. However, for the gluon transverse wave function we use the z-dependent gluon 
mass: mg{z < L) = rrig^p and mg{z > L) = rrig^y. 

Using the above formulas for the transition amplitude and wave functions with the 
help of the standard Fermi golden rule one can obtain after some simple calculations for 
the gluon distribution (below k = p^) 

dN rfiVW dN^^^ 



dxdli dxdh dxdk 



+ T3r2, (7) 



dxdk^ irx \ 2 y (k^ + /x2)2 



dxdk 



1-X + 



X 



TTX 



1 — COS 



(k' + f^l)L 



2E,x(l 



X] 



kV(fc)[(r(fc)-l)k' + r(fc)/ig-^g] 



(9) 
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where fx^ = rrigx'^ + — x) , r(k) = \Jas,p{k)/as^v{k). The first term on the r.h.s of 

(^) is the gluon spectrum in the vacuum, and the second one gives the medium correction 
of interest. In deriving (j^), (|^) we neglected the small spin-flip contribution (oc m"^). 

Recall that the gluon formation length is L/ ~ 2Eix{l — x)/(k^ + /i^) (we do not 
specify here the medium index since fip and /i^ are of the same order). Thus we see that 
the argument of cosine in the r.h.s. of @ is ~ L/Lf. Obviously, at -C L the rapidly 
oscillating cosine as function of L will vanish upon averaging over the production point 
of the fast quark, and one gets L-independent correction to the vacuum term. It can be 
written in the following physically transparent form 



dxdk 



dxdk 



+ ^^'^^ (10) 



where the first two terms describe simply modification of the spectrum due to the change 
in rrig and (in an infinite QGP), and the last term is the contribution of the transition 
radiation which reads 

dN^_CFa^/_ x^\ f kr(A;) k V 

dxdk^ ~ ^ V ^^^)W^u4 + j • ^ > 

It can be derived from (0) taking for the lower limit of the z-integral — oo. Note that the 
change in rrig and ag both cause the transition radiation. 

On the other hand, for the gluons with Lf ^ L expanding cosine in (^ one gets the 
correction oc 



dxdk^ 



Cpas{k)L' ( x'\ kMfc)[(r(A:)-l)k^ + r(fc)/ig-/i^ 



In this regime one cannot separate the transition radiation. Note that the above formulas 
demonstrate that the decrease of the cutoff and coupling constant in QGP work in opposite 
directions. One can also see that the relative contribution of the medium correction in (|^) 
is larger for gluons with Lf <L. These facts are consistent with intuitive expectation. 

4. Numerical results. In numerical calculations we take: m^ p = 0.4 and rrig y = 0.8 
GeV, for the quark mass we take m^^p = rng^^ = 0.3 GeV. As was mentioned above in 
the absence of accurate information on the in-medium as at low k we perform numerical 
calculations parametrazing it in the same F-model as for „ for several values of To 
understand the sensitivity of CFSI to a{p we use for it four values: 0.7, 0.5, 0.4, and 0.3. 
The first version corresponding to as,p = as,v is unlikely to be realistic since it neglects 
the in-medium modification of as{k), but it allows to see the magnitude of the pure mass 
effect. The last value is also unlikely to be realistic. Indeed, the results of the analysis of 
the lattice data within the quasiparticle model of QGP []T3[ say that as ~ 0.3 occurs for 



T ^ 2.2Tc. Since the QGP can have a temperature above this value only during a short 
time in the initial stage of its evolution the value of a{p = 0.3 is probably too small for 
evaluation of the CFSI. The value a(^„ = 0.5 seems to be most reasonable. For instance, 
it is close the value as ~ 0.47 obtained in ||15[ at T 1.5Tc ~ 250 MeV. Of course, it is 



not obvious that the thermal can be extrapolated safely to higher energies. However, 
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one can expect that in the sphtting of fast partons the effect of thermostat on as can only 
be weaker than that for the thermal partons. If this is the case, the CFSI correction may 
be larger than our estimate. 

In Fig. 1 we show the x- dependence of the ratio R{x) = ^ / '^^^."^ for several k- 
windows at Ei = 40 GeV averaged over L in the interval [0, 6] fm. As can be seen from 




Figure 1: The ratio R{x) = ^ / ^ at Ei = 40 GeV for a{; =0.7 (a), 0.5 (b), 0.4 
(c), and 0.3 (d) evaluated using Eqs. (7)-(9) for the /c-windows: [0,0.5] (solid line), [0.5,1] 
(dotted line), [1,1.5] (dashed line), and [1.5,2] (long dashed line) GeV. The ^-distributions 
obtained averaging over L in the interval [0,6] fm. 

Fig. 1 for tts^p = as,v the mass effect alone enhances considerably the probability of gluon 
emission at low x and k < 1 GeV. The curves for smaller values of a(^„ show that the 
enhancement of the radiation due to smaller cutoff in QGP is strongly compensated by 
the effect of smaller coupling constant in QGP, and for, probably unrealistic, a{p = 0.3 
there is a kinematic region where CFSI suppresses the gluon radiation. 
Using (H) we also calculated the energy loss, defined as 

AE = E, / dx / dk^x 



* dxd\<? 



For the limits of the x- and k^-integration we take Xmin = 'mg/Ei, Xmax = f^q/Ei, and 
^max = i[n.iia[Efx'^,Ef{l — x)^]. In Fig. 2 we show the results for AE as a function of L 
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at Ei =10, 20, 40, and 80 GeV. As one can see for af ^ =0.7, 0.5, and 0.4 the energy loss 




L, fm L, fm 



Figure 2: The quark energy loss AE due to the CFSI as a function of L for af ^ =0.7 
(a), 0.5 (b), 0.4 (c), and 0.3 (d) at Ei =10 (sohd line), 20 (dotted line), 40 (dashed line), 
and 80 (long dashed line) GeV. 

is positive and rising with L and Ei. For a{^p = 0.3 AE becomes negative which shows 
that suppression of the gluon radiation due to a small coupling constant becomes stronger 
than the enhancement caused by the mass effect. This strong cancellation between the 
two competing effects makes it difficult to give definitive predictions for the effect of CFSI 
in the region ~ 0.3 — 0.4. Note that our AE for af p = 0.5 appears to be of the 
same order of magnitude as the GLV prediction for the energy loss due to the induced 
radiation. 

5. Summary. We have shown that the final state interaction due to the change in 
the cutoff scale and running coupling constant from the vacuum to QGP modifies the 
gluon radiation from fast partons produced in AA-collisions. The contribution to the 
energy loss of this mechanism may be of the same order of magnitude as the induced 
gluon radiation. However, accurate evaluation of the CFSI is a difficult task since there 
are strong cancellations between the mass and running coupling constant effects, and the 
results depend strongly on the assumptions on the fc-dependence of the in-medium as{k). 

The results of the present paper raise a practical question, whether the jet tomographic 
analyses based on the theory of the induced gluon radiation can be used for extracting the 
density of hot QCD medium produced in AA collisions. At present, one cannot exclude 
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the possibility that the CFSI may appear to be even more important than the induced 
radiation. To clarify the situation it is highly desirable to study the influence of QGP on 
the running coupling constant for fast partons. 
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J. Speth for the hospitality at FZJ, Jiilich, where this work was completed. This work 
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